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In surface catalysis, the adsorption of carbon monoxide on transition-metal electrodes repre- 
sents the prototype of strong chemisorption. Notwithstanding significant changes in the molecular 
orbitals of adsorbed CO, spectroscopic experiments highlight a close correlation between the adsor- 
bate stretching frequency and equilibrium bond length for a wide range of adsorption geometries 
and substrate compositions. In this work, we study the origins of this correlation, commonly known 
as Badger's rule, by deconvoluting and examining contributions from the adsorption environment to 
the intramolecular potential using first-principles calculations. Noting that intramolecular anhar- 
monicity is preserved upon CO chemisorption, we show that Badger's rule for adsorbed CO can be 
expressed solely in terms of the tabulated Herzberg spectroscopic constants of isolated CO. More- 
over, although it had been previously established using finite-cluster models that Badger's rule is not 
affected by electrical conditions, we find here that Badger's rule breaks down when the electrified 
surface is represented as a periodic slab. Examining this breakdown in terms of anharmonic contri- 
butions from the effective surface charge reveals limitations of conventional finite-cluster models in 
describing electrical conditions at metal electrodes. 



I. INTRODUCTION 

The significance of carbon monoxide adsorption at 
transition-metal surfaces extends beyond its central rel- 
evance to catalytic exhaust control and electrochemi- 
cal energy conversion. Carbon monoxide is a corner- 
stone for fundamental research in surface science and 
electrochemistryi~— that has been sufficiently character- 
ized spectroscopically to serve now as a reliable ref- 
erence in measuring the potential of zero charge of 
metal electrodes^ - — identifying coadsorbates in com- 
plex surface environments,— ~— elucidating adsorption 
phenomena ; 11 ' 12 and probing electrical conditions in the 
electrochemical double layer£i2r— Among successful 
spectroscopic techniques, infrared experiments have de- 
livered quantitative insight into the properties of CO ad- 
sorbed at catalytic electrodes, unveiling distinctive vi- 
brational fingerprints as a function of the adsorption 
environment ,2iii"^2S In fact, the intramolecular frequency 
lj c of CO at close-packed metal surfaces undergoes red- 
shifts of several hundreds of cm' 1 with increasing ad- 
sorption coordination while augmenting gradually with 
the nobleness of the metal substrate^ In particular, this 
infrared trend has enabled for the vibrational recognition 
of CO adsorption sites on functional catalytic alloys^ 

Following the seminal work of Blyholder that has set 
forth the theory of electron donation (the depletion of 
the 5<r bonding frontier molecular orbital) and electron 
backdonation (the filling of the 2tt* and 2ir* antibond- 
ing orbitals) in discussing the electronic origins of CO 
chemisorption,— an important corpus of theoretical and 
experimental literature has developed on the elucidation 
of C-0 vibrational redshifts as a function of adsorption 
conditions. One salient, albeit unexpected spectroscopic 
feature that has attracted attention is the close correla- 
tion between the C-0 equilibrium stretching frequency 
w e and equilibrium bond length r c upon modifying the 



adsorption environment i i 25 In quantitative terms, the 
linear-regression correlation B — dui c /dr c (the Badger 
slope) is negative with an absolute value in the range 
6500-7500 cm _1 -A _1 regardless of the substrate compo- 
sition and adsorption site. In this work, we examine 
the electronic origins of this notable spectroscopic trend 
and carry the analysis further to probe the predictive ac- 
curacy of widely used first-principles computational ap- 
proaches in describing molecular adsorbates under elec- 
trical conditions. 

The study is organized as follows. First, we present 
structural and vibrational calculations for CO in the gas 
phase and focus on describing spectroscopic correlations 
for CO adsorbed on close-packed (111) metal surfaces 
from local and semilocal density-functional theories. Sec- 
ond, we introduce a quantitative intramolecular analysis 
to resolve anharmonic contributions from the surface en- 
vironment. On the basis of this analysis, we demonstrate 
that the correlation between oj c and r c for CO adsorbed 
at metal surfaces can be expressed in terms of the spec- 
troscopic constants of gas-phase CO. Finally, we explore 
implications of these results in simulating local electri- 
cal conditions and interpreting adsorption phenomena at 
electrified metal surfaces. 



II. FIRST-PRINCIPLES SPECTROSCOPY 

A. Carbon monoxide in the gas phase 

As a prelude to the intramolecular analysis presented 
in Sec. IIII1 we report electronic-structure predictions for 
the structural, vibrational, and nonlinear spectroscopic 
properties of CO in the gas phase. For the purpose of this 
analysis, particular attention is paid to the determination 
of high-order spectroscopic constants; those account for 
anharmonic contributions to the Born-Oppenheimer po- 
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TABLE I: Spectroscopic constants of CO compared with ex- 
periment (Ref. |26T). 





HF 


LDA 


PBE 


Expt. 


r c (A) 


1.087 


1.124 


1.141 


1.128 


cj c (cm x ) 


2434.14 


2128.59 


2103.93 


2169.81 


lo c x c (cm -1 ) 


11.47 


12.02 


11.91 


13.29 


Bo (cm -1 ) 


2.08115 


1.94558 


1.88862 


1.93128 


a e (cm -1 ) 


0.01597 


0.01774 


0.01778 


0.01750 


D c (lO^cm -1 ) 


6.0852 


6.5016 


6.08799 


6.1216 




FIG. 1: Comparison of intramolecular vibrational harmonic 
energy levels E v = hcw e (u + 1/2) and anharmonic energy 
levels E v = hc^2^L Y k0 (v + l/2) k for CO in the gas phase. 
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FIG. 2: Convergence of the Herzberg anharmonic spectro- 
copic constant a e for CO in the gas phase as a function of 
the order of the polynomial regression of the potential en- 
ergy curve (lower scale) and cutoff kinetic energy of the wave- 
function expansion (upper scale). 



tential energy V(r) (here, the coordinate r stands for the 
intramolecular bond length). 

In Fig. [TJ we compare C-0 vibrational energy levels 
evaluated within the harmonic approximation with the 
vibrational spectrum calculated taking into account an- 



harmonic contributions to the potential. Within the har- 
monic approximation, the potential energy can be ex- 
pressed in terms of the diatomic effective mass m, the 
equilibrium bond length r e , and the equilibrium stretch- 
ing frequency w e (in photon wavenumbers) as a quadratic 
function of the form 



V(r) = ^mc 2 (2TTu c ) 2 {r-r c ) 2 + 



(1) 



The diagonalization of the harmonic Hamiltonian yields 
vibrational energy levels that are separated by a constant 
energy gap, i.e., 



E„ 



hcw e (v + 1/2), 



(2) 



where v > denotes the vibrational quantum number. 

Considering now anharmonicity, the intramolecular 
potential can be written in the conventional adimensional 
Dunham form, 27 



V(0 = hca e(l + ai{; + a 2 e + ■■■), 



(3) 



where £ = (r — r c )/r c denotes the relative elongation 
of the bond and the terms ai are the coefficients of the 
Dunham expansion. Following the analytical treatment 
of Ref. [27], anharmonic energy levels can be evaluated as 



E„ 



+ OO 

hc}Y k0 {v 

k=0 



1/2)* 



(4) 



where the Yfco's are explicitly defined as a function of the 
Dunham et^'s. The first coefficient loo corresponds to the 
anharmonic energy shift while the coefficient Y10 can be 
identified straightforwardly to be 



Y10 — w c 



(5) 



The coefficient Y20 of the third term is typically negative 
and orders of magnitude smaller than Y10; it reflects the 
gradual narrowing of level separations in the upper part 
of the vibrational spectrum. Higher-order coefficients Y n Q 
are of alternating sign and decreasing magnitude. Quan- 
titatively, anharmonicity represents a negligible energy 
correction at the bottom of the vibrational spectrum, af- 
fecting the vibrational zero-point energy by less than a 
meV for gas-phase CO. In contrast, anharmonicity be- 
comes significant at higher vibrational energies. In fact, 
for v > 10, anharmonic corrections are on the order of 
one fourth of an eV, that is, comparable to the separation 
between energy levels. 

Alternatively to the Dunham adimensional represen- 
tation, in the Herzberg vibrational description, 28 anhar- 
monic levels are written in terms of the spectroscopic 
constants x c and y as 



E„ 



he 



+ 7^) _ u c x c {v + i) 2 
1, 



(6) 



3 



which allows us to identify Dunham coefficients as 



Y20 



-LO e X e 



(7) 
(8) 



Up to this point of the analysis, rotational degrees 
of freedom have not been considered. Quantized rota- 
tional energy levels labeled by the angular momentum 
quantum number J are obtained by adding the spherical- 
harmonics contribution 

J(J + 1) 



2mr 2 

to the vibrational Hamiltonian. The diagonalization of 
the rotationally augmented Hamiltionan yields the en- 
ergy levels 



+00 



E vJ = hcJ2 Y ki(v + V2)V(J + 1))', 



(9) 



k.l=0 



which are conventionally rewritten in the Herzberg rep- 
resentation as 



E u j = he 



+u cVc {v + i) 3 - a c (v + 1) J (J + 1) 
+B C J(J + 1) - D C (J(J + l)) 2 + ■ ■ 



(10) 



The coefficient a c that corresponds to the leading an- 
harmonic contribution and couples vibrational degrees 
of freedom to rotational degrees of freedom is of central 
significance to the analysis presented in Sec. IIIII 

To assess the performance of quantum-mechanical ap- 
proximations in describing the harmonic and anharmonic 
vibrational properties of CO in the gas phase, we perform 
first-principles spectroscopic calculations using the plane- 
wave CP (Car-Parinello) code of the QUANTUM-ESPRESSO 
distribution 29 that optimizes electronic degrees of free- 
dom via damped fictious Newtonian dynamics.— Our cal- 
culations are carried out at three different levels of quan- 
tum approximation, namely, uncorrelated Hartree-Fock 
(HF), the local density approximation (LDA)j^i and the 
semilocal Perdew-Burke-Ernzerhof (PBE) 32 approxima- 
tion with density-gradient corrections. The size of the 
cubic supercell is set to be of 30 bohr with countercharge 
corrections to eliminate dipole interactions between arti- 
ficial periodic images of the CO molecule.— In perform- 
ing gas-phase calculations, we employ norm-conserving 
pseudopotentials and select the cutoff energy of the wave- 
function plane-wave decomposition (smooth discretiza- 
tion grid) to be always higher than 30 Ry. Note that 
we use LDA pseudopotentials in our HF calculations, an 
approximation that has been justified by directly com- 
paring our HF structural and vibrational results to their 
tabulated counterparts^ In evaluating spectroscopic pa- 
rameters, we perform a polynomial regression of order 



8. With these calculation parameters, all of the spectro- 
scopic constants are found to be converged within less 
than a few percents. In particular, we verify that the 
predicted PBE anharmonic parameter a c is converged 
within 0.005 cm -1 , as shown on the two horizontal scales 
in Fig. H 

The results reported in Table U confirm the significance 
of DFT correlation in describing the structural and har- 
monic properties of CO. As a matter of fact, the pre- 
dicted LDA bond length r c — 1.124 A and vibrational 
frequency lu c = 2128.59 cm -1 are in closer agreement 
with experimental data, r e = 1.128 A and lu c — 2169.81 
cm -1 , than their HF counterparts, r c = 1.087 A and 
bj c = 2434.14 cm -1 (in accordance with the known HF 
propensity to overbind). The same trend is observed for 
anharmonic properties; including local electron correla- 
tion at the LDA level improves significantly the deter- 
mination of Herzberg constants, reducing the errors by 
nearly one order of precision in comparison with uncor- 
related HF. Semilocal PBE is also more precise than HF. 
However, PBE is found here to be slightly less accurate 
than LDA, exhibiting some tendency to underbind CO. 

Having confirmed the accuracy of gas-phase predic- 
tions in reproducing experimental data, we now proceed 
to assess the performance of LDA and PBE in describing 
adsorption geometries and vibrational modes for CO on 
select close-packed metal surfaces. 



B. Adsorbed carbon monoxide 
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FIG. 3: Equilibrium stretching frequency w e as a function 
of the bond length r c within LDA for CO adsorbed at the 
atop, bridge, and hollow sites on metal surfaces. The linear- 
regression slope is calculated to be -7659 cm _1 -A _1 . 

Herein, we focus on predicting structural and vibra- 
tional properties of CO adsorbed on close-packed sur- 
faces within LDA and PBE using the plane- wave PW code 
of the quantum-espresso distribution, which proceeds 
by iterative diagonalization of the effective Kohn-Sham 
Hamiltonian^ 

At this stage of the study, we emphasize that care must 
be exerted in applying LDA and PBE to describe CO 
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TABLE II: LDA and PBE structural parameters for CO adsorbed on close-packed (111) metal surfaces compared with experi- 
ment (r e denotes the equilibrium intramolecular bond length, r m stands for the adsorption distance from the carbon atom to 
the nearest metal atom, and 8 is the vertical tilt angle of the CO molecule). Unreported computational data correspond to 
bridge configurations undergoing a transition to the hollow site during structural optimization. 



Metal 


Site 


LDA 


PBE 


Expt. 


LDA 


'm {A) 

PBE 


Expt. 


a 
u 

LDA 


y ) 

PBE 


Ag 


atop 


1.140 


1.148 


— 


2.053 


2.173 


— 


0.8 


2.6 




bridge 


l.loo 


1. 101 




Z.LOO 


o 07c; 




n 7 
U. / 






hollow 


1.160 


1.167 


— 


2.230 


2.337 


— 


0.4 


0.3 


Au 


atop 


1.138 


1.148 


— 


1.960 


2.044 


— 


0.6 


2.4 






1 1 TO 


1 1 P.7 




2.082 


2.149 






9 s 
^ .0 




hollow 


1.169 


1.176 


— 


2.156 


2.228 


— 


0.4 


0.4 


Cu 


atop 


1.145 


1.154 


— 


1.822 


1.874 


1.91(1)" 


1.3 


1.2 




Krin (TO 
Ul JUtL 




1 1 7D 






2.002 






n r i 

VI . 




hollow 


1.168 


1.178 


— 


2.000 


2.062 


— 


0.7 


0.6 


Pd 


atop 


1.143 


1.152 


— 


1.849 


1.891 


— 


0.4 


0.2 




bridge 




















hollow 


1.174 


1.184 


1.15(5)'' 


2.040 


2.078 


2.05(4)' 


0.2 


0.1 


Ft 


atop 


1.144 


1.154 


1.15(5) c 


1.837 


1.865 


1.85(10) c 


0.4 


0.7 




bridge 


1.167 


1.177 


1.15(5) c 


2.001 


2.032 


2.08(7) c 


1.2 


1.8 




hollow 


1.178 


1.188 




2.084 


2.120 




0.4 


0.4 


Rh 


atop 


1.148 


1.159 


1.20(5)"* 


1.825 


1.855 


1.87(4) d 


0.7 


0.3 




bridge 


1.169 


1.180 


1.15(10) e 


1.988 


2.022 


2.03(7) e 


1.1 


1.3 




hollow 


1.177 


1.188 


L^IO/ 


2.059 


2.097 


1.99(7/ 


0.3 


0.0 



"ARPEFS (angle-resolved photoemission extended fine structure technique, Ref. [35]). 
6 LEED (low-energy electron diffraction technique, Ref. l36h . 
TEED (Refs.|2iand[37|). 
d LEED (Ref.il). 

TEED (bridge structure obtained in the presence of near-top coadsorbed CO, Ref. |39TI. 
TjEED (hollow structure obtained in the presence of coadsorbed benzene, Ref. I40T) . 



cal DFT approximations occasionally fail in predicting 
the relative stability of adsorption sites. These quali- 
tative deficiencies take root into spurious electron self- 
interaction (i.e., the incorrect analytical behavior of the 
approximate ground-state energy as a function of the 
electron number, which favors adsorption sites of high 
coordination). One notorious instance of the LDA and 
PBE failures is the overstabilization of multifold adsorp- 
tion sites (hollow and bridge) relative to terminal ad- 
sorption sites (atop) for CO at close-packed platinum 
surfaces4£~— 

In spite of quantitative errors in predicting the adsorp- 
tion energetics, it has been established from the analysis 
of the density of states (DOS) and force density of states 
(FDOS) that semilocal DFT approximations are accu- 
rate in determining structural and vibrational properties 
for adsorbed COi22 The calculations presented below cor- 
roborate the performance of LDA and PBE in describing 
structural and vibrational properties. 

Our computational benchmark consists of a pool of 6 
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FIG. 4: Equilibrium stretching frequency ai e as a function 
of the bond length r c within PBE for CO adsorbed at the 
atop, bridge, and hollow sites on metal surfaces. The linear- 
regression slope is calculated to be -7151 cm _1 -A _1 . 



adsorption. In fact, it is known that local and semilo- 
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TABLE III: LDA and PBE intramolecular vibrational fre- 
quencies lj c (cm -1 ) for CO adsorbed on close-packed (111) 
metal surfaces compared with low-coverage experimental 



data. 


Metal 


Site 


LDA 


PBE 


Expt. 


Ag 


atop 


2084 


2026 


2042(6) a 




bridge 


1965 


1913 


1970(10)° 




hollow 


1919 


1875 


1884(21)° 


An 


atop 


2107 


2035 


2060° 




bridge 


1938 


1882 






hollow 


1859 


1817 


— 


Cu 


atop 


2073 


2009 


2075 c 




bridge 




1873 






hollow 


1887 


1825 


<1840 c 


Pd 


atop 


2103 


2033 


<2090 d 




bridge 










hollow 


1851 


1786 


1840 d 


Ft 


atop 


2118 


2048 


2070 e 




bridge 


1917 


1851 


1830 e 




hollow 


1825 


1759 


1760 e 


Rh 


atop 


2079 


2004 


1990^ 




bridge 


1899 


1832 


<1870 / 




hollow 


1835 


1771 





"FTIR (Fourier transform infrared spectroscopy, Ref . l4ll ) . 
"PR-IRAS (polarization-modulation infrared 
reflection-absorption spectroscopy, Ref. y). 
TRAS (infrared reflect ion- absorption spectroscopy, Ref. I42T). 
d SFG (sum- frequency generation spectroscopy, Ref. l43h. 
e SFG (Ref. 03). 

'HREELS (high-resolution electron energy loss spectroscopy, 
Ref.H). 



representative metal elements with face-centered crystal 
structures, namely, Ag, Au, Cu, Pd, Pt, and Rh. First, 
we determine lattice parameters and bulk moduli at the 
different levels of DFT approximation, controlling con- 
vergence with respect to the cutoff kinetic energy of the 
plane-wave expansion of the electronic wave functions 
(and total electronic density when employing ultrasoft 
potentials), to the Monkhorst-Pack sampling of the Bril- 
louin zone, and to the temperature of the generalized 
electronic entropy- 
Predicted crystal properties are reported and com- 
pared with experiment in Appendix |XJ Based upon cal- 
culated lattice parameters, we construct periodic (111) 
surface models adopting a fully relaxed orthorhombic 
(-\/3 x 2) four-layer slab geometry to determine the vi- 
brational properties of CO adsorbed at the atop, bridge, 
hollow (face-centered cubic) sites, corresponding to a 
low coverage of 1/4 ML. We underscore that adopt- 
ing a model x 2) geometry is not restrictive for 



the present analysis as the resulting variations in struc- 
tural and vibrational predictions are negligible on the 
scale of the Badger correlation as discussed further be- 
low. All of the considered adsorption systems corre- 
spond to upright, nondissociative adsorption and are 
not affected by spin-polarization. In determining and 
analyzing vibrational frequencies, we employ a frozen- 
phonon method restricted to vertical C and O atomic 
displacements. In this frozen-phonon picture, one deter- 
mines the stretching frequency from the largest eigen- 
value of the two-dimensional intramolecular dynamical 
matrix Dij — (m,imj)~^ d Zi d Zj V(z±, Z2) where Z\ and z^ 
denote the vertical displacements of the C atom of mass 
mi and O atom of mass m^- The restriction to trans- 
verse adsorbate coordinates represents a reliable approx- 
imation for the intramolecular vibrational mode of light 
CO on heavy metals, yielding equilibrium stretching fre- 
quencies in close agreement with the result of the full 
DFPT (density-functional perturbation theory)^ calcu- 
lation of the density matrix, with a maximal error of 2% 
for CO sitting at the most coordinated (i.e., hollow) site 
of the lightest (i.e., Cu) substrate. 

In Table UH we report LDA and PBE structural pa- 
rameters, namely, the equilibrium bond length r c , the 
distance from the carbon atom to the closest metal atom 
r m , and the angular tilt of the CO molecule 6. Our 
computational predictions are compared with available 
experimental data, mainly derived from low-energy elec- 
tron diffraction (LEED). Of note in Table ILTl is the lack 
of experimental structural data for CO adsorbed on no- 
ble metals and the difficulty to probe adsorption sites 
of high coordination without specific coadsorption [as is 
the case for CO at the hollow site of Rh(lll)]4& Notwith- 
standing experimental uncertainties and variations in ad- 
sorption configurations, LDA and PBE predictions al- 
ways lie within (or very close to) the experimental range. 
Precisely, LDA tends to underestimate r e by a few hun- 
dredths of A, whereas PBE exhibits improved predictive 
performance. Additionally, LDA and PBE calculations 
confirm the experimentally observed upright adsorption 
geometry for all of the surfaces considered. 

In Table [TTTT local and semilocal intramolecular vibra- 
tional frequencies are presented and compared with low- 
coverage experiments based upon infrared techniques, 
including infrared reflection-absorption spectroscopy 
(IRAS), sum- frequency generation (SFG), and high- 
resolution electron energy loss spectroscopy (HREELS). 
LDA and PBE calculations predict important vibrational 
rcdshifts as a function of the adsorption coordination, as 
large as of 300 and 345 cm -1 for the hollow site of Pt. In- 
stead, the lowest redshifts are observed for CO adsorbed 
at the terminal site on Pt. From available experimen- 
tal frequencies, LDA and PBE mean absolute errors are 
estimated to be of 40 and 30 cm -1 , respectively. 

With these results in hand, we can graphically confront 
predicted C-0 vibrational frequencies u> c with equilib- 
rium bond lengths r e , thereby recovering the close corre- 
lation highlighted in Ref. I25I . In quantitative terms, lin- 



6 



ear regression yields Badger slopes of —7659 and —7151 
cm _1 -A _1 within LDA and PBE, respectively (Fig.[3]and 
Fig. HI). It is important to note that our PBE prediction 
is in very good agreement with the slope of —7066 cm -1 
reported in Ref. [2l| using the semilocal PW91 (Perdew- 
Wang) approximation^ with a different adlayer model, 
confirming the validity of our structural approximation. 
In the next section, we present a quantitative procedure 
to rationalize Badger's rule in terms of gas-phase spec- 
troscopic properties. 

III. ANALYSIS 
A. Analytical procedure 



V (eV) 




FIG. 5: Anharmonic stretching direction (dotted line) in 
the plane of mass-scaled atomic coordinates (y'miZi, y/rn^Z2) 
for CO at the atop adsorption site on Ag. The anharmonic 
stretching direction connecting energy minima departs from 
its harmonic counterpart (dashed line with arrow). 

In a nutshell, the intramolecular analysis developed 
here consists of deconvoluting the intramolecular poten- 
tial energy along the curved anharmonic stretching tra- 
jectory of CO in adsorption (corresponding to the dotted 
line in Fig. [3]). 

To determine the anharmonic stretching trajectory, we 
first calculate the potential energy surface correspond- 
ing to constrained vertical displacements of the C and O 
atoms, the total energy being minimized with respect to 
all of the other molecular and metal-atom coordinates. 
For instance, we consider the case of CO adsorbed at 
the atop site of Ag in Fig. [5] The two-dimensional po- 
tential energy surface is sampled with 17 points unifor- 
mally distributed along the Zi, z 2 , {i\ + z 2 )/a/2 and 
(zi — z 2 )/-\/2 unit vectors that define the C, O, C-O, 
and Pt-CO displacement directions. As already men- 
tioned, for light CO adsorbed on heavy metals, restrict- 



ing the analysis to transverse adsorbate coordinates pro- 
vides an accurate description of the vibrational behavior 
of CO with errors that do not exceed 1-2% in predicting 
the intramolecular stretching frequency. The stretching 
trajectory is obtained by connecting the points derived 
from the minimization of the total energy fixing the in- 
tramolecular distance. Minimum-energy points are cal- 
culated by simple steepest descent, projecting the energy 
gradient on the direction defined by the constraint on 
the bond length. Finally, the total energy is plotted as a 
function of the constrained bond length, thereby obtain- 
ing the intramolecular potential energy curve describing 
the elongation of adsorbed CO. 

This analytical construction provides a straightforward 
generalization of harmonic modes to the anharmonic case 
and offers a direct means to extract the intrinsic vibra- 
tional properties of any polyatomic system coupled to 
a substrate. Note that the curved stretching trajectory 
reduces to a straight line when the potential energy sur- 
face is represented by a parabolic well, in agreement with 
physical intuition. As a consequence, the curvilinear de- 
viation of the stretching trajectory reflects anharmonic 
contributions to the intramolecular potential energy of 
relevance to the present analysis. 



B. Intramolecular correlation 
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FIG. 6: Intramolecular potential Vkds(r) for CO adsorbed at 
the atop site of Pt, intramolecular potential V ga s(j") in the gas 
phase, and their difference V a ds(f) — V sas (r). 

Using the procedure outlined above, we now examine 
the contribution from the adsorption environment to the 
intramolecular potential. For this analysis, we focus on 
semilocal PBE, which has been shown above to be more 
accurate than local LDA in predicting the vibrational 
properties of adsorbed CO. 

In Fig. HI the intramolecular potential V a d s {r) for CO 
at the atop site of Pt is compared to the intramolecular 
potential V ga , s (r) for CO in the gas phase. Evaluating 
the difference between the two potentials, we obtain the 
adsorption contribution to the intramolecular potential, 
which is close to linear on a large range spanning all of 
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FIG. 7: Badger slope B^ds (left scale) and percent deviation 
of the Badger slope (right scale) relative to its gas-phase coun- 
terpart as a function of the bond length r c for CO adsorbed 
at the atop, bridge, and hollow sites of metal surfaces. 



the bond lengths r e reported in Table |TTJ This provides 
computational evidence for the preservation of the anhar- 
monicity of the intramolecular potential despite strong 
electronic hybridization. 

To extend this observation, we now consider the other 
adsorption systems described in Sec. Ill Bl In Fig. [7] (left 
scale), we calculate the Badger slope B a d s at the mini- 
mum of the intramolecular potential: 



B a ds = 



ads 



dr 



(11) 



where uj a d s (r) is related to the second derivative of the 
potential through 



' ads 



dr 2 



mc 2 (271-Wads) 5 



(12) 



We first observe that the influence of the substrate drops 
as the activity of the metal and coordination of the ad- 
sorption site decrease; the contribution from the sub- 
strate is found to be maximal at the hollow site of Rh 
and Pt while is minimal at the atop site of Ag and Au. 
Despite these fluctuations, the magnitude of the slope 
£> a ds remains within 6400-7600 cm -1 -!" 1 , explaining the 
correlation between w e and r c depicted in Fig. 2) 

To gain further quantitative insight, we compare cal- 
culated Badger slopes B a d s to the Badger slope of the 
isolated molecule 



B ri 



du„ 



dr 



(13) 



where ^J gas (r) is related to the second derivative of the 
potential in the gas phase. The slope £> gas can be calcu- 
lated from the potential energy V gas (r), yielding a cor- 
relation coefficient of —7601 cm _1 -A _1 . Alternatively, 
B gaB can be evaluated by exploiting the Dunham expan- 
sion [Eq. (JT}]. As a matter of fact, making use of the 
relation 



(l + oi) + 



(14) 



that is valid to leading order in the adimensional num- 
ber B 2 /uj 2 [see Eq. (15) in Ref. 2J7j, one arrives at the 
analytical expression 



B„ 



3ujc 

2r~. 



1 



~6BT 



(15) 



Substituting the terms in Eq. (1151) for the calculated val- 
ues reported in Table [3J we obtain a slope of —7600 
cm _1 -A _1 , which is in very good agreement with the 
value of —7601 cm _1 -A _1 determined above. By com- 
paring this slope to that of adsorbed CO [Fig. [7] (right 
scale)], we conclude that the magnitude of the surface 
contribution does not exceed 16% and is even lower than 
5% for terminal sites on noble metals. 

The central conclusion of this analysis is that the hy- 
bridization of CO molecular orbitals with metal bands 
causes a maximal departure of 1200 cm _1 -A _1 from the 
gas-phase Badger correlation. This trend provides quan- 
titative justification of Badger's rule and allows us to 
express the slope B a ds for CO in adsorption in terms of 
the anharmonic constants of the isolated molecule. In 
the next section, we re-examine this conclusion taking 
into account the influence of electrical conditions and we 
discuss the validity of conventional adsorption models in 
describing electrified surfaces. 



IV. DISCUSSION 
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FIG. 8: Intramolecular frequency cj c as a function of the 
bond length r c for CO adsorbed at the atop, bridge, and 
hollow sites of Pt considering a varying surface charge. De- 
parture from Badger's rule is expressed in terms of percent 
deviations of the charge-dependent slope B c \ cc relative to the 
gas-phase slope Z3 gas . The dashed line serves as a guide for 
the eye. 

The intramolecular analysis has demonstrated the pos- 
sibility of embedding the complexity of the adsorption en- 
vironment into a simple linear constraint that preserves 
Badger's rule. Based upon this observation, it is tempt- 
ing to infer that Badger's rule is also conserved under 
electrical conditions. Here, we confront commonly used 
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adsorption models, namely, finite-cluster and periodic- 
slab simulations under electrical conditions to discuss the 
validity of this conjecture. 

We first concentrate on finite-cluster models, which 
have been employed in simulating adsorption in the low- 
coverage limit i 25 ' 50 " — To analyze the influence of elec- 
trical conditions within finite-cluster models, we exploit 
data from Ref. [HI, which have been computed at the 
scmilocal generalized-gradient level using a finite Pti3 
cluster consisting of 2 hexagonal layers of 7 and 6 atoms; 
such a Pti3-cluster structure allows us to study atop ad- 
sorption on the 7-atom facet and hollow adsorption on 
the 6-atom facet. In these calculations, the influence of 
the electric field is accounted for via a transverse lin- 
ear potential added to the effective electronic Hamilto- 
nian and interatomic potential. Based upon extracted 
data (cf. Fig. 2 and 3 of Ref. l5lh , we calculate the field- 
dependent slope 



dwc 
~dE 



E=0 



dr c 
dE 



E=0 



(16) 



to be 6790 cm" 1 ^" 1 at the atop site and 6180 cm" 1 ^" 1 
at the hollow site. The calculated slopes are in very good 
agreement with the slopes £? a ds reported in Sec. MI Bl 
(Fig. [7]). In quantitative terms, B c \ cc deviate by less 
than 5% from B a & a for both atop and hollow adsorptions. 
These results provide evidence that applying a uniform 
electric field within finite-cluster simulations conserves 
Badger's rule. 

Having verified that Badger's rule is preserved within 
finite-cluster models, we now turn to periodic-slab 
calculations M^— In periodic-slab calculations, electri- 
cal conditions can be imposed by applying an exter- 
nal electric field E similarly to finite-cluster calculations 
(the electric-field method) or, alternatively, by varying 
the surface charge a through adjusting the number of 
electrons in the system (the surface-charge method). 62 
For periodic slabs of sufficient thickness, the two ap- 
proaches can be shown to be equivalent >§2 Here, we adopt 
the surface-charge approach with density-countercharge 
techniques 59 to impose two-dimensional boundary condi- 
tions compatible with the slab periodicity. We consider 
slabs comprising 3 metal layers and verify that increas- 
ing the thickness to 4 layers affects Badger-slope pre- 
dictions by less than 150 cm" 1 '!" 1 . The magnitude of 
the surface charge is raised up to 0.1 C-m" 2 i££ As de- 
picted in Fig. [51 we calculate charge-dependent slopes 
^cicc by simultaneously monitoring the vibrational fre- 
quency ui c and the intramolecular distance r c . We thus 
obtain charge-dependent slopes of —5700 cm _1 -A _1 at 
the atop site, —5070 cm" 1 -A" 1 at the bridge site, and 
—5010 cm" 1 !" 1 at the hollow site. Consequently, cal- 
culated slopes are found to deviate significantly from 
their gas-phase counterparts, demonstrating that Bad- 
ger's rule does not hold under electrical conditions using 
periodic-slab models at variance with finite-cluster sim- 
ulations. 



To explain these trends, we resort to effective electri- 
cal variables, which provide a conceptually useful, direct 
mapping for describing the influence of electrical condi- 
tions on physical systems. Explicitly, we first write the 
free energy of the adsorption system as 



U(r,E) = V(r) + SP{r,E)E 



(17) 



in terms of the polarization of the adlayer per unit sur- 
face area P = dU / d(SE). Variations with respect to the 
intramolecular distance r yields an implicit electric-field 
dependence for the equilibrium bond length of the form 
dV/dr\ r=rD{E) + Ser*(r c (E), E)E = where 



<9P 

dr 



(18) 



is known as the Born effective charge per surface area 
and is commonly used in describing the coupling between 
structural properties and an external electric field. 48 
Making use of this equilibrium relation, we derive the 
expression of the field-dependent Badger slope 



B 



clcc 



Buds — 



lo c 9 In | a* 
2 



dr 



(19) 



t— r c ,E— 



which shows that B c \ cc departs from S a ds by a term that 
is related to the logarithmic derivative of the effective 
charge <J*(r, E) as a function the bond length. 

As a result, the clear accordance between B e i ec and 
B a ds within cluster models indicates that the effective 
charge of the adlayer remains constant upon stretching 
the CO adsorbate. In other words, within finite-cluster 
descriptions, the dipole P(r, E) depends linearly on the 
bond length r due to the absence of charge redistribution 
at the surface of the cluster. 

In notable contrast to finite clusters, a nonlinear cou- 
pling mechanism takes place on periodic slabs; substitut- 
ing our PBE predictions into Eq. (fl9l) . we can assess the 
coefficient of the leading term 7 that describes the redis- 
tribution of the adlayer charge upon elongating the CO 
molecule, 



1+7- 



(20) 



where a* — a*(r Cl E = 0). The adimensional coefficient 
7 is calculated to lie between —2 at the hollow site and 
— 1.5 at the atop site, reflecting the nonlinear dependence 
of the surface dipole as a function of the bond length. 
This is at variance with finite-cluster models that fail to 
capture the decay of the effective charge. In conclusion, 
finite-cluster models are not apt at describing the struc- 
tural and vibrational properties of CO adsorbates un- 
der varying electric field, strongly suggesting that they 
should be employed with caution in simulating electrical 
conditions at charged electrodes in contrast to periodic- 
surface models that correctly describe the decaying effec- 
tive surface charge of the adlayer. 
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V. CONCLUSION 

In summary, Badger's rule, which correlates linearly 
the equilibrium stretching frequency uj c to the bond 
length r c of CO under varying adsorption conditions 
arises from the conservation of intramolecular anhar- 
monicity despite strong chemisorption. The intramolec- 
ular analysis has shown that the hybridization of CO 
molecular orbitals with the metal electronic states causes 
moderate deviations from the gas-phase Badger correla- 
tion for all of the metal electrodes considered; in partic- 
ular, Badger deviations become negligible for adsorption 
at the terminal site of noble metals. 

Applying a uniform electric field within finite-cluster 
models has been shown to not affect Badger's rule sig- 
nificantly, whereas the anharmonic contribution from an 
electrical conditions within periodic-slab models has been 
found to alter the intramolecular correlation, ultimately 
leading to the breakdown of Badger's rule. Analyzing 
the influence of electrical conditions in terms of effective 
electrical variables reveals that finite-cluster electric-field 
models are not apt at capturing variations in the effec- 
tive surface charge upon elongating the CO adsorbate, 
explaining their failure in recovering the influence of elec- 
trical conditions on the intramolecular bond. 

This first-principles study underscores the relevance 
of the redistributed surface charge in understanding the 
coupling of the adsorbate vibrational modes with the 
surface electric field and points out limitations of finite- 
cluster models in simulating adsorption under electrical 
conditions. These results motivate further theoretical 
and computational efforts towards more realistic simula- 
tions of electrode interfaces — ideally taking into account 
the self-consistent polarization and ionic response of the 
embedding electrolyte^ Such computational approaches 
would find fruitful applications in studying molecular 
electrochemical phenomena relevant to catalysis and en- 
ergy conversion. 
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Appendix A: Crystal data 

For completeness, we report calculated lattice param- 
eters and bulk moduli together with experimental data 
in Table IIV1 We observe that lattice parameters are un- 
derestimated by LDA with an absolute error of approxi- 
mately 2% and overestimated by PBE with an error that 
does not exceed 3%. Despite uncertainties in experimen- 
tal bulk moduli, errors in predicting bulk moduli are ad- 
mittedly nonnegligible, reaching 35% for Ag within LDA 

TABLE IV: LDA and PBE lattice parameters and bulk mod- 
uli for Ag, Au, Cu, Pd, Pt, and Rh in the face-centered-cubic 
crystal structure compared with experiment (Ref. |26T). 



DFT Expt. 



Metal 




a (A) 


B (GPa) 


a (A) 


B (GPa) 


Ag 


LDA 


4.019 


135.9 


4.086 


100 




PBE 


4.163 


90.7 


Au 


LDA 
PBE 


4.048 
4.173 


191.1 
140.0 


4.078 


200L20 


Cu 


LDA 


3.552 


171.6 


3.615 


140 




PBE 


3.674 


127.7 


Pd 


LDA 
PBE 


3.864 
3.969 


224.7 
169.7 


3.890 


180 


Pt 


LDA 


3.921 


298.6 


3.924 


230 




PBE 


3.995 


246.8 


Rh 


LDA 
PBE 


3.784 
3.859 


310.5 
259.3 


3.803 


327.5L52.5 



and 30% for Au within PBE. Although inaccuracies in 
predicting the elastic properties of the substrate may af- 
fect the description of vibrational modes involving dis- 
placements of metal atoms, the CO stretching mode is 
weakly coupled to the substrate and is thus moderately 
affected by such errors, as shown in Sec. Ill Bl 
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A wave-function cutoff energy E cut = 30 Ry is sufficient 
to achieve targeted convergence levels for structural, har- 
monic, and anharmonic properties with the exception of 
Cu (40 Ry for LDA and 35 Ry for PBE), Pd (40 Ry for 
PBE), Pt (40 Ry for PBE), and Rh (35 Ry for LDA). In 



our ultrasoft calculations, the charge-density cutoff energy 
is set to be 8_E C ut. For bulk crystals, the Monkhorst-Pack 
sampling of the Brillouin zone is set to be 10 x 10 x 10 (se- 
lecting a Marzari- Vanderbilt generalized-entropy tempera- 
ture of 0.03 Ry) with the exception of Au (LDA) and Cu 
(PBE) for which the Brillouin resolution has to be raised 
to 12 x 12 x 12. Based upon calculated lattice parameters, 
an equivalent surface density of sampling points is applied 
in adsorption calculations. 

1 C-m~ 2 of surface charge corresponds to 11.29 V-A -1 of 
vacuum electric field. 



